Bifurcation delay - the case of the sequence: 
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^ Abstract 

f~] ' Let us give a two dimensional family of real vector fields. We suppose that there exists a stationary point 

. where the linearized vector field has successively a stable focus, an unstable focus and an unstable node. 

When the parameter moves slowly, a bifurcation delay appears due to the Hopf bifurcation. The studied 
0^ i question in this article is the continuation of the delay after the focus-node bifurcation. 
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1 Introduction 



"Singular perturbations" is a studied domain from many years ago. Since 1980, many contributions were written 
because new tools were applied to the subject. The main studied objects are the slow fast vector fields also 
known as systems with two time-scales. We will give the problem here with a more particular point of view: 
the bifurcation delay , as in articles [8, 2, 9, 7]. We write the studied system: eX = f{t,X,e), where £ is a real 
positive parameter which tends to zero. For a better understanding of the expression dynamic bifurcation it is 
I better to write the system after a rescaling of the variable: 



X = f{a,X,e) 

a = e 



where a is a "slowly varying" parameter. 

The main objects in this study are the eigenvalues of the linear part of equation X = f(a,X,0) near the 
^*)' \ quasi-stationary point. Indeed, they give a characterization of the stability of the equilibrium of the fast vector 
field at this point. The aim of this study is to understand what happens when the stability of a quasi-stationary 
point changes. A bifurcation occurs when at least one of the eigenvalues has a null real part. 

In this article we restrict our study to two-dimensional real systems. In this situation, the generic bifurcations 



H \ are: the saddle-node bifurcation, the Hopf bifurcation and the focus-node bifurcation. 

The saddle-node bifurcation is solved by the turning point theory: when the real part of one of the eigenvalue 
becomes positive, there is no delay and a trajectory of the systems leaves the neighborhood of the quasi- 
stationary point when it reaches the bifurcation. For this study, the study of one-dimensional systems is 
sufficient: we have a decomposition of the phase space where only the one-dimensional factor is interesting. 
There exist many articles on this subject, we will be interested particularly by [3] where the method of relief \s 
used. The article [6] introduces the geometrical methods of Fenichel's manifold. 

The Hopf delayed bifurcation is well explained in [10], we will upgrade the results in paragraph 2 below. 
In a focus-node bifurcation, the stability of the quasi stationary point does not change, then, locally, there 
is no problem of canards or bifurcation delays. Indeed, when there is a bifurcation delay at a Hopf-bifurcation 
point, it is possible to evaluate the value of the delay, and the main question is to understand the influence of 
the focus node bifurcation to this delay. 

In paragraph 2, the Hopf bifurcation alone is studied, as well as the focus-node bifurcation following a Hopf 
bifurcation in paragraph 3. 
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In paragraphs 2.1 and 3.1, we assume that there exists a solution of the system approximed by the quasi 
steady state in the whole domain, so this trajectory has an infinite delay. The used methods are real, and the 
system has to be smooth (actually only C^). In paragraphs 2.2 et 3.2, we avoid this very special hypothesis. It 
is here supposed that the system is analytic, and we study the solutions on complex domains. Unfortunately, I 
have not a proof for the main result of this article. But it seems to me that the problem is interesting, and the 
results are argumented. 

We use Nelson's nonstandard terminology (see for example [5]). Indeed, almost all sentences can be translated 
in classical terms, where e is considered as a variable and not as a parameter. Often, the translation is given 
on footnotes. 

2 The delayed Hopf bifurcation 

The problem is studied and essentially resolved in [10]. We give here the proofs to improve the results and to 
fix the ideas for the main paragraph of the article. The main tool is the relief's theory of J.L. Callot, explained 
in [4]. 

The studied equation is 

eX = f{t,X,s) (1) 
where / is analytic on a domain X» of C x x (D. 

Hypothesis and notations 

HI The function / is analytic. It takes real values when the arguments are real. 
H2 The parameter e is real, positive, infinitesimal^. 

H3 There exists an analytic function 4>, defined on a complex domain Vf so that f(t,(p{t),0) = 0. The curve 
X = (pit) is called the slow curve of equation (1). We assume that the intersection of "Dj with the real 
axis is an interval \tm,tM[- 

H4 Let us denote \{t) and ^j.{t) for the eigenvalues of the jacobian matrix Dxf , computed at point {t, (f){t), 0). 
We assume that , for t real, the signs of the real and imaginary parts are given by the table below : 



t 




a 


tM 


SR(A(i)) 









^{^i{t)) 









9(A(i)) 












h + 





Then, when t increases from tm to tM, the quasi-steady state is first an attractive focus, then a repulsive 
focus, with a Hopf bifurcation at t = a. 

2.1 Input-output function when there exists a big canard 

In this section, we assume that there exists a big canardX(t) i.e. a solution of equation (1) such that^ ^{t) ~ ^(<) 
for all t in the S-interior of ]tm, tni- We now want to study the others solutions of equation (1) by comparison 
with X. 

The main tool for that is a sequence of change of unknowm: first, we perform a translation on X, depending 
on t to put the big canard on the axis: 

X = X{t) + Y 

It gives the system 

sY = g{t, Y, e) with g{t, Y, s) = f{t, X{t) +Y,e)- f{t, X{t), e) 

^In classical terms, we assume that e leaves in a small complex sector: |£| bounded and arg(£) e] — S,S[. 

^Without nonstandard terminology, a big canard is a solution of equation (1) depending on the parameter e such that 

Vt6]tm,tM[, lim X{t,e) = <l>{t) 

£>0,£—^0 
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The matrix Dx.f(t. 4>(t).0) has two complex conjugate distinct eigenvahies (see hypothesis H4), then there exists 
a Hnear transformation P{t) which transforms the jacobian matrix in a canonical form. We define the change 
of unknown 

Y = P{t)Z 

The new system, has the following form (we wrote only the interesting terms): 



eZ = h{t,Z,e) with h{t,Z,e) 



a{t) -uj{t) 



Lo{t) a{t) 

The next change is given by the polar coordinates: 



Z + 0{e)Z + 0{Z^) 



X{t) = a{t) - ioj{t) 



rcos9 
r sin^ 



sr = r {a{t) + 0{e) + 0{r)) 
e0 = u{t) + 0{e) + 0{r) 

The last one is an exponential microscope^: 



r = exp 



p = a{t) +0{e) +e^ki{r,0,e) 
eB = Lj{t) + 0{e) + eJk2{r,9.e) 



(2) 



While p is non positive and non infinitesimal, r is exponentially small and the equation (2) gives a good 
approximation of p with p = a. When p becomes infinitesimal, with a more subtle argument (see [1]) using 
differential inequations, we can prove that r becomes non infinitesimal. This gives the proposition below: 

Proposition 1 Let us assume hypothesis HI to H4 (Hopf bifurcation) for equation (1). However, we assume 
that there exists a canard X{t) going along'^ the slow curve at least on \tm,tM{- Then if X{t) goes along the 
slow curve exactl'if' on]te,ts[ with [te,ts] C]tm,tM[, then 



i: 



SR(A(T))dT = 



The input-output relation (between t^ and tg) is defined by J*" iR(A(r))<iT 
(see figure 1). In this case, this relation is a function. 



0. It is described by its graph 




Figure 1: The input-output relation for equation (3) when there exists a big canard. 



^All the preceeding transformations were regular with respect to e. This last one is singular at e = 0. 
*A solution X{t,e) goes along the slow curve at least on ]ti,t2[ if 

Vte]ti,t2[, ^lim ^ X(t, e) = </)(<) 

solution X{t,e') goes along the slow curve exactly on ]te,ts[ if it goes along the slow curve at least on ]te,ts[, and if the 
interval ]te,ts[ is maximal for this property. 
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2.2 The bump and the anti-bump 

In this paragraph, t becomes complex, in the domain Vt. We assume that for all t in Vt, the two eigenvalues 
X(t) and are distinct. It is a necessary condition to apply Callot's theory of reliefs. 
We define the reliefs Rx and i?^ by: 




Rx{t) = ^{Fx{t)) 



J a 

It is easy to see that A(t) = and F\{t) = F^{t), then Rx(t) = Rfj.(t). The two functions Rx and R^ 
coincide on the real axis. We will denote R{t). 

Definition 1 We say that a path 7:56 [0, 1] i-^ X>t goes down the relief Rx if and only if -j^Rxilis)) < for 
all s in [0, 1] . 

Definition 2 Let us give a point t^, such that (ie,0(ie),O) G V. We say that Vt is a domain below if and 
only if for all t in the S -interior ofVt, there exist two paths in Vt, from te to t, the first one goes down the 
relief Rx and the second one down R^. 

Theorem 2 (Callot) Let us assume that Vt is a domain below t^. A solution X{t) of equation (1) with 
an initial condition X{te) infinitesimally close to (j^ite) is defined at least on the S-interior of Vt where it is 
infinitesimally close to 4>{t). 

Let us apply this theorem to the following example, chosen as the typical example satisfying hypothesis HI 
to H4 (Hopf bifurcation) . 

f ei ^ tx + y + eci . . 

\ ey = -X + ty + ec2 

The eigenvalues are A = t—i et /i = t+i. The level curves of the two reliefs Rx{t) = and Rp,{t) = ^{t+i)'^ 

are drawn on figure 2. 




Figure 2: The level curves of the two reliefs of equation (3), and a domain below ie 

Generically^ there is no surstability at point t = i (see [3] for the definition of surstability) . Consequently, 
we have the following results for all equations such that the reliefs are on the same type as those of figure 2: 

Definition 3 Let us give tc a point^ where the eigenvalue vanishes: A(tc) = 0. The value of the relief at point 
tc is a critical value of the relief Rx- The bump^ is the real number t* bigger than a, minimal such that Rx{t*) 
is a critical value. The anti-bump is the real number t** smaller than a, maximal such that Rx{t**) is a critical 
value. 

do not know the exact generic hypothesis. We have to combine the constraints given by the surstability theory of [3] and the 
fact that the equation (1) is real 

^In some cases, it is possible that tc is infinite. For eX = ^ ^^ost sint ^ ^ ^ 0{X^) + 0(e) we have tc = +ico. 
^The name "bump" is a translation of the french name "butee" 
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For equation (3), the bump is = 1 and the anti-bump t** = —1. 

Theorem 3 A trajectory of equation (1) can go along the slow curve X = (f){t) exactly on \te,ts[ if and only if 
one of the following is verified: 

te < t** and ts = t* 
te = t** and ts > t* 
t** <te < a and a < ts < t* and R{te) = R{ts) 

This theorem is illustrated by the graph of the input-output relation, drawn on figure 3. 



te 




Figure 3: The input-output relation for equation (3) 



3 Delayed Hopf bifurcation followed by a focus-node bifurcation 

The studied equation is 

eX = f{t,X,e) (4) 
where / is analytic on a domain D of (D x (D^ x (D, and satisfies the following hypothesis: 

Hypothesis and notations 

HFNl The analytic function / takes real values when the arguments are real. 
HFN2 The parameter e is real, positive, infinitesimal. 

HFN3 There exists an analytic function <±>, defined on a complex domain Vf such that f{t,(j){t),0) = 0. The 
curve X = (j)it) is called the slow curve of equation 1. We assume that the intersection of Vt with the real 
axis is an interval ]tm,iM[- 

HFN4 Let us denote X{t) and ii{t) for the eigenvalues of the jacobian matrix Dxf , computed at point (t, <j){t),Q). 
We assume that , for t real, the signs of the real and imaginary parts are given by the table below : 



t 




a 


b 




U{X{t)) 







h + 


+ 


W)) 







h + 


+ 


3(A(t)) 














+ 


+ 


h 






Then, when t increases on the real interval ]tm,tM[, we have succesively an attractive focus, a Hopf 
bifurcation at t = a, a repulsive focus, a focus-node bifurcation at t = b and a repulsive node. At point 
t = b, the two eigenvalues coincide. We assume that X{t) = ii{t) only at point b. Actually, in the complex 
plane, the two eigenvalues are the two determinations of a multiform function defined on a Riemann 
surface with a square root singularity at point b. 

However, there is a symmetry: if the function ^ is defined with a cut-off on the positive real axis, it 
satisfies \fs = — \/i and we then have 

ii{t) = Mt) 
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HFN5 For the same reason, the two rehefs 

Rx{t)=^(^J X{t)dt^ and Ri^{t) = ^(^J n{t)d?j 

are the two determinations of a multiform function with a square root singularity at point t = b. However, 
there is a symmetry: if the function ^ is defined with a cut-off on the positive real axis, it satisfies 

^/s = — \/5 and we have then: Rix{t) = R\(t) except on the cut-off half line [6, -|-oo[. For real t > b, 
we choose determinations of square root such that X{t) < iJ,{t). We assume that R\ has a unique critical 
point with critical value Re- We assume that Rx{b) < Re- An example is given and studied in paragraph 
3.2.1. 



3.1 Input-output function when there exists a big canard 

We assume now that there exists a big canard X{t) i.e. a solution of equation (1) such that X{t) ~ (j){t) for 
all t in the S'-interior of ]tm,tM[- The study below is similar to paragraph 2.1. The added difficulty is the 
coincidence of the two eigenvalues at point b which do not allow to diagonalize the linear part. 
The first change of unknown is X = X{t) + Z which moves the big canard on the axis X = 0: 

sZ = A{t)Z + 0{e)Z + 0{Z^) , A{t) = Dxf{t,(l>{t),0) 
Let us denote ^ ^ the coefficients of the matrix A{t). As in paragraph 2.1, the change of unknowns 



rcos9 \ ( p\ 

rsme ] ' ^ = ^"PUJ 



gives the new system: 



p = a{t) cos'^0 + (Pit) + -fit)) cose smd + d{t) 81-0^9 + 0(e) +eeki{r, 9, s) 
eO = Jit) cos"^ + {6{t) - a{t)) COS0 sm9 - I3{t)sm^ 9 + 0{e)+esk2{r,9,s) 



(5) 



= arctan 



For nonpositive p (more precisely, for infinitesimal r), the second equation is a slow-fast equation. Its slow curve 
is given by 

'S{t) - a{t) ± ^J a{tf - 2 a{t)5{t) + + 4 /3(t)7(i) \ 

m) ) 

It has two branches when A and /x are reals, one is attractive, the other is repulsive: see figure 4. 

When 9 goes along a branch of the slow curve, (and when r is infinitesimal), an easy computation shows 
that p is infinitely close to one of the eigenvalues A or p. The repulsive branch corresponds to the smallest 
eigenvalue (which is real positive). When t < b, the angle moves infinitely fast, and an averaging procedure 
is needed to evaluate the variation of p: 

i:r''id0 

An easy computation shows now that, in the S-interior of the domain t < b, p < 0, we have 

ip) - ^^^1±M = 3?(A(i)) = 5i(/,(i)) 

Let us give an initial condition (t, 9) between the two branches of the slow curve and p negative non 
infinitesimal (in the example, we can take t = 0.8, = 0, p = —0.03). For increasing t, the curve {t,9(t)) goes 
along the attractive branch of the slow curve, while p believes negative non infinitesimal. For decreasing t, the 
solution goes along the repulsive branch, then moves infinitely fast while p believes negative non infinitesimal. 
Consequently, we know the variation of p{t) (see figure 4). As in paragraph 2.1, a more subtle argument is 
needed to prove that when p becomes infinitesimal, the variable r becomes non infinitesimal and the trajectory 
X leaves the neighborhood of the slow curve. 

From this study, all the behaviours of p{t) are known, depending on the initial condition. They are drawn 
on figure 5. 
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Figure 4: One of the trajectories of system 0.002X 



t 1 
t-0.3 t 



) 



X drawn with the variables {0,p). The 



slow curve is also drawn 




Figure 5: The possible behaviours of p{t). 



Proposition 4 Let us give an equation of type (6) with hypothesis HFNl to HFN5. Assume also that there 
exists a big canard X{t) going along the slow curve on the whole interval ]tm,tM[- If a, trajectory X{t) goes 
along the slow curve exactly on an interval ]te,ts[ with [te,ts] c]tm,^A/[, then 



Conversely, if the inequalities above are satisfied, there exists a trajectory going along the slow curve exactly 

on ]te,ts[. 

The input-output relation is described by its graph, drawn on figure 6. 

We could give more precise results if we consider the two variables r and 6 for the input-output relation. 
Indeed, when the point {te,ts) is in the interior of the graph of the input-output relation, we know that, at time 
of output, 9 is going along the attractive slow curve which corresponds to the unique fast trajectory tangent to 
the eigenspace of the biggest eienvalue ^. 

3.2 The focus-node bifurcation is a bump 

Here is the main part of this article. Today, I am not able to prove the expecting results, but I have propositions 
in this direction. To explain the problem, I will give conjectures. 

Let us define the anti-bump t** and the two bumps t*^ and i* as in definition 3: 




Rx{t**) 



Rx{t\) 



R^^it;) 
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Figure 6: The input-output relation for equation (6) when there exists a big canard. 

. We have t** < a < t*^ — t*^ < b or t** < a < b < t*^ < t\. In the first case, the bump is before the 
focus node bifurcation, and the study of paragraph 2.2 is available. The interesting case is the second, where 
the computed bump is after the focus node bifurcation, this case is assumed with hypothesis HFN5. 

Conjecture 5 With hypothesis HFNl to HFN5, the following proposition is generically wrong: 

If a trajectory of (4) goes along the slow curve at least on ]t**,a[, then it goes until the slow curve at least 

on 

To work on this conjecture, we will study an example which is, in some sense, a normal form of the problem: 
the slow curve is moved on the i-axis and the fast vector field is linearized. The example is 

r = tx + y + e^ci , , 

\ e^y = {t- b)x + ty + e^C2 

Proposition 6 A numerical simulation of equation (6) gave the figure 7. It confirms conjecture 5. 





1 


0.2 

-X 

b 


t 1 




1 ^—-^ 





Figure 7: Trajectories X_ and X^: the first goes along the horizontal axis from — oo to b, where it jumps outside 
the neighborhood of the horizontal axis; the second one goes along the horizontal axis from +00 to —b where it 
has big oscillations. The parameters are b — 0.3, ci — 0, 02 — —1, = 0.002, the trajectory is computed with 
a RK4 method, with step 0.0001. Other methods and other steps were tried, and the results are always very 
similar. 

This proposition gives a good argument for the next conjecture, more precise than the first one: 

Conjecture 7 If a trajectory of system (4) goes along the slow curve in a neighborhood of a real t with t < a 
and R{t) > R{b), then it does not go along the slow curve after the focus-node bifurcation point b. 
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So, generically, the input-output relation of equation (4) has a graph similar to the graph of figure 3; if R{f**) > 
R{b), we have to replace t* by b et t** par t^* where R{tl*) = R{b). The delay of the Hopf bifurcation is stopped 
either by the bump (as in case of a Hopf bifurcation alone) either by the focus-node bifurcation. 

Proposition 8 // the conjecture 7 is true for one trajectory, then it is true for all of them. 

Proof Assume that equation (4) has a solution X which does not verify conjecture 7. Then, X goes along 
the slow curve on an interval ]ti,f2[ with t\ < tl* < a < b < t2- If the problem is considered on a restricted 
interval ]ti,t2[, the equation has a big canard, and we can apply the proposition 4. Then ah trajectories going 
along the slow curve before t^* goes along the slow curve until b, and even a little more. ■ 

In this article, we will now study only equation (6). We changed e into only to avoid fractionnary 
exponents. The analytic structure with respect to e is obviously modified, but does not matter for our purpose. 

To study the phase portrait of equation (4) or (6), two trajectories are very important. They are called 
distinguished trajectories by JL.Callot and they are very classical. The first one, denoted X+ goes along the 
slow curve for t near tM- Similarly, X_ goes along the slow curve for t near tm- These two trajectories are 
Fenichel's manifolds, they are unique when tm = —oo and Im = +oo. For the particular equation (6), these two 
trajectories are drawn on figure 7. We have for this example a nice fact: X_ and X+ have an explicit formula, 
using the Airy function (in an appendix (section 4) , we give classical needed results on Airy functions and Airy 
equation) . 

X+(i) = ) = -e'^^Mit)£'^e-i^M-\r)dr(^ ^ ) (7) 

where M{t) = ^ ( sj^A'l^-^) ) ^^^^^^^^^ = I 

All the integrals are convergent because the Airy function is bounded at infinity by C|t|~^e3l*l^ . 
3.2.1 The relief 

In this paragraph, wc want to explore the methods used in paragraph 2.2 when there is a focus-node bifurcation. 
We also check the hypothesis HFNl to HFN5. 

Hypothesis HFNl to HFN3 are obvious with the slow curve 4){t,X,0) = and the domain D = (D x x C. 

The computation of the eigenvalues of the jacobian matrix J{t) = ^ ^ ^ ^ | ^ gives 

X(t) = t - {t - b)i iJ,{t) = t + {t - b)i 

The determination of the square root is needed to ahow the formula above. In all this paragraph, we choose 
a cut-off on the positive real axis: 

{re'")^ = V? 0G [0,2tt[ 

For the function we choose the same cut-off. 

The relation = —t^ will be useful. Then, A and fj, are the two determinations of a multiform function. 
The cut-off is the semi-axis [b, +oo[, and /u(t) = X(t). 
For a = and 

b>i, (10) 

the hypothesis HFN4 is easy to check. 
The two associated reliefs are given by 

Fx{t) = ht" - i(* - b)i - libi F^{t) = ii^ + f (i - b)i + iibi 

Rx{t)=n{F^{t)) R^{t)=di{F^{t)) 
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Figure 8: Level curves of relief R\ for b = 0.3, and path used in paragraph 3.2.4. 



Let us comment figure 8: the value of Rx is +oo at both ends of the real axis. If a path goes from t — — oo 
to t = +00, it has to go down at least until the mountain pass, which is the unique critical point of the relief 

given by 

= 1 +iy^b- i = 0.500 + 0.224 i 

The value of the relief at this critical point is 

Rc^Rxitc) = ib-J^ = 0.067 

We solve now on the real axis the equation Rx{t) = Rc. The solution are and ts given by 

if 6>i + iV3 , ie = ~Vb^ ts^^/b^ 
if 6<i + iv^ , te = -Vb^ { Z - 

The symbols ... in the formula above are the solutions of a polynom in t of degree 4. The exact expression is 
not needed. For 6 = 0.3, we have 



te 



-0.365 



tsi = 0.346 ts2 = 0.525 



The value tsi is on the sheet right to the cut-off: arg(tsi) = 2tt. Besides, ts2 is on the sheet left to the cut-off: 
arg{ts2) = 0. When we look on the polynom which has tsi and is2 as roots, we can prove that the hypothesis 
HFN5 is satisfied for 

1 < b < i + iVs (11) 



3.2.2 Callot's domains 

To study the canards of equation (6), we introduce two special solutions, called distinguished solutions by J.L. 
Callot: = has an asymptotic^ condition X^{+oo) = and X_ — has an asymptotic 

condition cx)) = 0. They are unique. In this paragraph we build a domain I?+ where X+ is infinitesimal 

(it corresponds in the complex plane to the expression "going along a real interval"). In allmost all situations, 
the builded domain is the maximal domain with this property. 

^Here the things are easier than in the general case because the domain Vt contains the whole real axis. In general case, there 
is no unicity of the distinguished solution, but the difference remains exponentially smaller than the computed quantities. 
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For trajectory X+ In this paragraph, it is better to change the cut-off, and we define (only in this paragraph) 

We are looking for a complex domain 1)+ such that the real point +00 is in I?+ , the singularity b is not in I?+ . 
We look for domains below +00 (see definition 2) for the relief R\ and also below +00 for the relief R^. 

On figure 9, such domain is drawn^ in dark. Attention: at the left, the domain has a spike with a real part 
smaller than —b and a nonzero imaginay part. The intersection of I?+ with the real axis is ] — &, 6[U]6, -|-c»[. 
The theorem of Callot (theorem 2) says that is infinitesimal on the whole S'-interior of 

Actually, a more precise study shows that the domain is not the maximal domain where X+ is infinites- 
imal: if we consider domains on the the Riemann surface (two sheets covering) we can add to I?+ its conjugate 
(drawn in lightgray on the figure 9). Because the solution X+ is analytic without singularity at point 6, it is 
infinitesimal on the symetric domain. 





■ ■f' ' , • ' 
■^■v.' V"^'. "-' .- . 

>■",-■ .-■ 


/ / 1 i / » = i 5" -I 

— v/ ■i:E-- 

0T^_ Hi Q.« I "D-GO ■ 11.75 B D.ai 


- - J\ i 


! : - I ■ ■ 

\ i "i = \ 




y \ li 



Figure 9: The domain I?+ for 6 = 0.3 



For trajectory X- A similar method gives the domain V- such that X- is infinitesimal on the S'-interior of 
T>- . It is easier because we do not need to consider a two sheets covering. The domain I?_ is drawn on figure 
10. 

3.2.3 Evaluation of X+(&) 

The slow curve a; = y = is repulsive for all positive t. Then the trajectory X+ is infinitesimal at least for all t 
positive non infinitesimal (in fact it is infinitesimal on a larger interval). Its asymptotic expansion in power of 
is given by formal identification in the equation: X^ — X]n>o ^n(^)^^" has to verify the recurrence identities 

Vn-\ = (^ - &)a;„ + <y« + '5„-iC2 

where 8n~\ = 1 if n = 1 and vanishes for all others n. The computation of the first terms is easy: 

_ -cii + C2 3 t(^^+t-3b)ci + (-3t2 + t + 6)c2 ^6 9 

_ (t - &)ci - tea 3 (-2^3 -H 35^2 + &t - b2)c^ + (^3 + -m-t + b)c2 g 9^ 
and we we have now proved the 



^The picture is a little bit different when b is greater or smaller than — f — i\/123. For this particular value, we have R^i{b) 
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Figure 10: The domain I?_ for 6 = 0.3 



Proposition 9 



X+{b) 



( 



) 



3.2.4 Evaluation of X (6) 

The simple method above is not convenient to evaluate X-{h) because we expect that X- does not go along 
the slow manifold in a neighborhood of b. 

We will use the explicit formula (8) to evaluate X_(6). The computation is a little bit tedious. In all the 
formulae below, the symbol represent a quantity which goes to zero when e > goes to zero. 

The inverse of the matrix M is easy to compute: we know the determinant of M (see the property 4 in the 
appendix on Airy's functions). 



To compute the integrals in formula (8), we change the real path of integration ] — oo, b]. For some integrals 
we choose a path which goes down the relief R\ from — oo to b, for other integrals, we choose the conjugate 
path which goes down the relief (the idea is the same as in Callot's proof of theorem 2). The path which 
goes down Rx is drawn on figure 8. The end of the path is a vertical segment from 6 + i/3 to b. At point 6, it is 
tangent to the level curve of the relief, then, the path does not go down the relief with the precise definition 1. 
Thus, we have to be care with approximations at this point. 




Let us denote 




It is one of the function we have to integrate to evaluate X- . 



Lemma 10 Let us give r such that t — b is non infinitesimal and ^tt < arg(T — 6) < tt. Then 




Proof Using the asymptotic expansion of A (see in appendix), we have: 




) 



4 



(l + 0(e^)) 
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Substituting in the formula of /, we have: 

eHn\f{T)\ = SR(i62-iT^-f(/(T-6))i+0) 

We write t — 6 in polar coordinates: t — b = re*^, with 9 G]|7r, 7r[. Then j^(T — b) = re^^^^^^K Because 6' — §7r 
has an argument between — jTT and ^tt, the power | gives (j^(T — 6)) ^ = g^i^o-^) _ This expression can be 
writed — (t — with the same determination of as in i?A- ■ 

The interesting consequence of this lemma is that along the considered path, the function / is increasing with 
a logarithmic derivative of type s~^. To precise, we need the following lemma: 

Lemma 11 There exist two constants k and S standard? , positive such that 



p. ^ - ^ 

£2 



V<7e[0,^] , \f{b + iae^)\<ke-^''' 



Proof By definition of /, we have 



f{b + iae^) = e-'^'''ei'''''A{ifa) then \f{b + iae^)\ = e^"^ " \A{if 
For real positive infinitely large a, the asymptotic expansion of the Airy function give the estimation 



2^ 



e 3^ 



(we know that 5R((y^)2) = Then if is real standard, less than we have the following inequality, 

true for all a infinitely large: 

\A{ifa)\ < e-'^'^^ 

By permanence", this inequality beheves true for all real a greater than some positive standard iv. We can 
deduce the following majoration: 

Vae[a;,^] , \f{b + iae'')\ < e^'^^"^^'^^ 

For a < CO, we have: 

V£7e[0,w] , \f{b + iae'^)\ < ei'^^^ki < 2ki with fci = max \A{ifa)\ 

iT£[0,a>] 

Then we are looking for a constant S such that 

VcT G [0, ^] , iCT^e - (5icr5 < -Sai 

The inequality is equivalent to cr < K^L^i^ \ choice of 6 less than 5i — is convenient. This choice is 

possible only if Si > |\//? what is true as soon as /? < | and 5i near enough from 
To verify the majoration of the lemma for a < co, we can choose 

k = 2fcie^'^' 



The next lemma is the more technical part of the article. The purpose is to evaluate an oscillating integral 
with successive integrations by parts. 



■'here, it is tiie same to assume tliat k and 5 are independent of e 

■^The non standard arguments in these proofs can be replaced by classical arguments, but, for that, new quantified variables 
have to be added, and it seems to me that the idea of the proof is more understandable with nonstandard language. 
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Lemma 12 We have the following expansion: 



/ 

Jb 



Proof Let us substitute r by 6 + ie^a in the integral. We have 

[ f{T)dT = -ie^ 1^ f{b + ie^a)da = -ie'^ e-i'''ei''^'A{ifa)da 

Jb+il3 Jo Jo 

The exponential e"^'^' is fast oscillating. The exponential e^"^ ^ is infinitely close to 1 for all non infinitely large 
a and A{ij'^(j) is decreasing. All properties are checked to apply the method of integrations by parts. But there 

1 2 

is a difficulty: e^"^ ^ is increasing and does not believe close to 1. Now, let us explain the computations: 

I = [ fiT)dT = h+h + h with 

h = j{f{b + m-f{b)) h = - J Jj^ af{b + ie^a)da 

/3 = f{b + is'^a)da with f{b + is'^a) = e-^^'e^"^' A'{ifa) 

Jo 

With lemma 11, we know that f{b + (ii) is exponentially smaller than f{b) = A(0). Thus, we have 

7i = _lA(O)e3 + 0£7 
To estimate I2, we perform a new integration by parts: 

h = Ji + J2 + J3 + J4 with 



ie'' (3 ie' 



72- 



= ~^^fi^ + f^') = "-yr a^f{b + ie''a)da 

J3 = Jj^ af{b + ie^a)da J4 = '-^ Jj^ f{b + ie^a)da 

Because f{b + Pi) is exponentially small, we have Ji = 0e^. We have also J4 = — fr-f- If you substitute A' for 

.2 

A the expression I3 is the same as All the arguments are the same with function A' and function A. Let 

-V .s -2 ^ 

us denote /j, Ji the expressions obtained from li and Ji when ^' is substituted for A. Thus we have I3 = 
To estimate J2 we perform a new integration by parts exactly as for evaluation of I2 ■ J2 = Ki + K2 + K3 + K4. 
All the integrals are bounded by a non infinitely large real number because all the integrated functions are 
bounded (see lemma 11) by a integrable standard function. To summarize: 

I = h+h + h l2 = Ji + J2 + J3 + Ji J2 = Ki + K2 + K:i + Ki 

7i = _iA(O)£3 + 0e^ Ji = pe"^ Kr = 

l2 = ^e^ J2 = K2 = -^ J ' a^f{b + ie^a)da = 

/3 = Y^ J3 = ^72 i^3 = ^J2 -^4 = -^/ Ki = —^l2 (12) 

Then, all the ingredients are given, and we can compute the asymptotic expansion of / in powers of e. To start, 
we have I — (J)s. For similar reason, I = (/)s. Then, using formulae 12, we have I3 = then / = 0£^. We 
iterate the process, inserting the known approximations in formulae 12, and we obtain a better approximation: 
Is = (pe^ then 

I = -Ia{0)s^ + 
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The next step: 



The next step: 

/3 = -^A'(O)e4 + 0£4 J3 = 0e4 j^ = ^^(o)e'^ + ^£6 l^ = (f)e^ 

I = -^A(0)£3 _^A'(0)e4 + 
The next step: (do not use the relation ^"(0) = 0, because we have sometimes to substitute A' for A): 

h = -^^A'{O)s^-^A"{O)e^ + 0s^ J3 = l2 = 0e^ 

I = -^A(0)e3 -^^A'{0)e^ - ^A"{0)e'' + 

/3 = -^^'(O)£^-p^"(O)£'-^A"'(O)£6 + 0e6 J3^0£6 i;C3 = £'^ 

J4 = lA(0)£6 + iJA'(0)£^ + (Z)£7 if4 = (Z>£^ 
J2 = 0£6 72 = ^^(0)£6 + p£6 

7 = -iA(0)£3 -■^A'iO)e^ - ^A"iO)s' + (^1a(O)-1a"'(O))£'^ + 0£« 

The last step: 

73 = -^^A'{0)e''-^A"iO)e'-^A"'{0)e'+(^^A'{0) ^^""(0)) £^+0£^ J3 = ^A'(O)£^+0£^ 7^3 = 
7 = - Ja(0)£3 - gA'(0)£4 - ^A"{0)s' + (1a(0) - ^A'"(0)) £'^ + (^A'(O) - ^A""(0)) £^ + 0£^ 



Lemma 13 



J f{T)dT= e^p(^-^iRx{b + ip)-Rx{b) + 0)^ 

Proof The chosen path goes down the relief 7?;^, then the lemma is a corollary of the majoration of lemma 
10. ■ 

Lemma 14 

l\i'^A'(j^^-^yr = -iA'(0)£3 - ^^(0)£^ + (1-^)a'(0)£« + 
f_J^'^A'{j"-^)dr = -Ja'(0)£3 - g^(0)£^ + (^_^)^'(0)£e + 0£^ 
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Proof With lemmas 12 and 13, the firs part of the lemma is proved. A similar computation giATS the second 
part, if we remember that A"(0) = but A'"{0) ^ 0. So, the vanishing terms are not the same in the two 
formulas. The two last formulas are the complex conjugate of the two first one. ■ 



Proposition 15 



X-{b) = 



Proof Insert the estimations of lemma 14 in the explicit formula (8), and, after tedious simplifications, the 
proposition is proved. ■ 

Conjecture 16 The two values X-{b) and X+{h) have the same asymptotic expansion. 

With Maple, I checked that the two expansions coincide until terms in . 



4 Appendix: Airy's functions 

The Airy's equation is linear, non autonomous of second order. It is 

^=tx (13) 

The pair {A{t), B{t)) of Airy's functions is a fondamental system of solutions. The function satisfy the following 
properties (these results can be found in every book on special functions). 

1. The value at the origin are: 

^(0)=3-i-i- A'(0) = - — ^ B(0) = 3-5-i- B'iO) = —^ 

^ ' r(|) ^ ' 2 TT ^ ' r(|) ^ ' 2 TT 

2. On a sector of angle less than §7r, around the positive real axis^, the Airy's functions have an asymptotic 
expansion for t going to infinity: 

Ait) = J-e-i*^-i(l + 0(t-i)) A!{t) =-J_e-i*^ti(l + 0(t-i)) 

B{t) = ^ei*^t-i(l + 0(t-i)) B'{t) = ^ei*^ti(l + O(i-i)) 
The functions A ei B are oscillating when t goes to — oo. 

3. Let us denote j = e^^^ = — 3 + -^i- The Airy's equation is invariant by the change of variable t jt, 
then A{jt) and B{jt) are also solutions. So they can be written as a Hnear combination of A{t) and B{t). 
We perform an identification at point to find the coefficients: 

A{jt) =-ifA{t) + ^ifBit) B{jt) = iifA{t)-\fB{t) 

Mft) = -hiA{t)-li3B{t) B{jH) - -iijA{t)-^jB{t) 

4. Classicaly, the couple {A{t),B{t)) is chosen for a base of the set of solutions. It could be better (in a 
study in the complex plane) to choose {A{jt), A{j'^t)) for base. With Liouville's theorem, we prove that 
the following determinant is constant, and we compute its value at the origin. 



det ' ^(^'^^ ^(■^■'*) ^ - A 



ft _ 
jA'ijt) fA'ift) J - 2n 



^Take care: the determination of is here the classical determination with a cut off on the negative real axis, not the 
determination choose along all this article 



16 



Figure 11: Graphs of Airy's real functions A and B 
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